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Abstract. A numerical study of quenched QCD for light quarks is pre- 
sented using 0(a) improved fermions. Particular attention is paid to the 
possible existence and determination of quenched chiral logarithms. A 'safe' 
region to use for chiral extrapolations appears to be at and above the 
strange quark mass. 



1. Introduction 

The goal of lattice QCD is the computation of physical quantities such as 
hadron masses and matrix elements using numerical Monte Carlo methods. 
This has proved to be an ambitious programme because after discretisation 
of the path integral and generation of a sufficiently large number of inde- 
pendent configurations several limits must be considered: 
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1. The box size. This is currently at ~ 1.5 — 3 fm, and should be compared 
with the nucleon rms radius of ~ 0.8 fm. 

2. The chiral limit m q — > 0. The u/d and s quarks are light quarks. 

3. Continuum limit a k — > (k = 2 if we choose O(a) Symanzik improved 
fermions, staggered fermions or Ginsparg- Wilson fermions; for Wilson 
fermions we expect the discretisation effects to have k = 1). 

If all these limits can be successfully taken then presumably QCD will 
reproduce nature. Although first attempts in this direction are being made, 
[1], it will require much faster computers to achieve this goal. To reduce the 
computational effort often the quenched approximation is employed when 
the fermion determinant is simply set to a constant. However then new 
problems arise (or are exacerbated): 

1. Spurious quenched chiral logarithms appear as m q — > 0. 

2. The appearance of exceptional configurations. 

3. Consistency of the final results when comparing with their experimen- 
tal or phenomenological values. 

In this talk we shall consider points 1 and 2 numerically using 0(a) im- 
proved fermions. 

2. Chiral perturbation theory and quenched chiral logarithms 

This was developed by Bernard and Golterman, [2] and Sharpe, [3]. What 
do we expect? The quenched pseudoscalar effective chiral Lagrangian gives 

{amps) oc (am q ) i+« , 

a 2 g P = (0\V\ps) oc (K) 2 ] 4 , 

am ps af ps = (Q\Ai\ps) oc am ps . (1) 

As the rj remains light and has a single and double pole in its propagator, 
this latter term at p 2 = 0, (wio/ m ps)(V m ps) ac ^ s n ^e an extra vertex giving 
a singular correction to the usually harmless loop term m ps ln(m ps /A 2 ) of 
Trig ln(m ps /A) 2 , so that the logarithmic term becomes singular, [4]. This 
can then be summed to give eq. (1). Normally PCAC, d-A = 2fh q V, would 
give us (am ps ) 2 oc afh q . Thus the non-zero 5 leads to singular behaviour 
for m ps and gp in the chiral limit for quenched QCD. Simple estimates 
lead to an expectation for 5 ~ 0.1 — 0.2. A is a cut-off on the r/' loop so 
A < 900 MeV and above this scale any singularities are surely damped out. 
For consistency, we would expect little difference between the PCAC quark 
mass, afhq, and the 'standard' quark mass, am q = |(1/k — l/^ c ). 

Similarly by considering the vector/baryon effective chiral Lagrangians, 
it can be shown that [5, 6] 

am v ,N = Co + c\am ps + c 2 (om ps ) 2 + 0((am ps ) 3 ), (2) 
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where there is an extra linear term present, as the r/ gives to the usual 
term (am ps ) 3 an additional term m^am ps . 

3. Numerical results 

How do these theoretical considerations fare with the numerical data? Nu- 
merically, it advantageous to use the PCAC quark mass, fh q because this 
quark mass can be found very accurately and does not depend on the 
first-to-be-determined parameter k c . It is convenient to consider the ratio 
afhq/ (am ps ) 2 . This is shown in Fig. 1, [7], for degenerate quark masses at 
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Figure 1. arh q /(am ps ) 2 against (am ps ) 2 for 0(a) improved fermions. Filled circles 
denote points used in the linear fits. The dashed line is the mass of an (unphysical) 575s 
pseudoscalar meson, (using tuk), while the dotted line represents m x . 

values of 6.0, 6.2 and 6.4. To give an idea of scales, we note that using the r$ 
'force scale' then m, lies almost at the chiral limit (within our numerically 
accuracy there is no difference between these points) and a hypothetical 
pseudoscalar meson, composed of the strange quark and its antiquark, lies 
at about (rom ps ) 2 ~ 3.13 (when using uik)- For the charm quark (using 
ttid) we find (rom ps ) 2 ~ 44.9, way off the plot scale. 

It is to be seen from the picture that from the strange quark mass 
to heavier quark masses we have linear behaviour. (Indeed the linearity 
seems to hold until rather heavy quark masses, say m q < \m c .) Below the 
strange quark mass, there seems to be a (sharp?) break in this behaviour. 
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Possibly we can attribute this to the onset of quenched chiral logarithms. 
mp S = \/2mK corresponds here to about 700 MeV ~ A so one might hope 
that any quenching effects are suppressed above this value. We now check 
the behaviour between afh q and am q . In Fig. 2 we plot am q /am q against 
am q . While there seems to be a reasonably linear relation between the 
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Figure 2. am q /arn q against am q forO(a) improved fermions. Filled circles denote points 
used in the fits. 

two (lattice) definitions of the quark mass above am s , below we again see 
deviations. This fit is somewhat sensitive to the value of k c used; although 
the general picture shown in Fig. 2 never seems to change significantly. 
(Indeed using all the light quark data in the fit still produces a similar 
result.) This perhaps obscures the interpretation of Fig. 1 as being due to 
quenched chiral logarithms. Nevetheless, due to problems in determining 
am q , we prefer to use the results with am q , [8]. 

To try to expose the small quark region, and to determine S (if the 
deviations are due to quenched chiral logarithms) then it is convenient to 
plot the logarithm of Fig. 1. This is done in Fig. 3, where we expect 
the slope to be 5 (for the fitted quarks masses below the strange quark 
mass). We find for [3 = 6.0, 5 ~ 0.12(4), and for (3 = 6.2, 5 ~ 0.06(2) (for 
(3 = 6.4 there is not enough data). For (3 = 6.0, at least, there is reasonable 
agreement with the theoretical prejudice; for (3 = 6.2 the value seems small. 

Despite the above results it should be noted, as discussed above, that 
it is notoriously difficult to numerically detect quenched chiral logarithms. 
Indeed the above effects may simply be due to finite-size effects or 'excep- 
tional configuration' problems. We have only been able to check very few 
quark mass points for finite size effects. The impression is that they are 
small; this is backed up by [8], who work on a larger lattice. 

Exceptional configurations are seen as either the non-convergence of the 
fermion matrix inversion or the correlator seems to have a 'fake source' at 
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Figure 3. \n[(afh q /(amn) 2 ] against ln[((ro/a) 2 x (arn ps ) 2 ]. The left line is a linear fit to 
the quarks with mass below the strange quark mass, while the dotted line is the previous 
fit from Fig. 1. 



some t value. The problem is more severe for 0(a) improved fermions than 
for Wilson fermions and increases as (3 [ and/or c sw | and/or m q j. (This 
is the main obstacle for the 0(a) improvement programme going below 
about (3 ~ 6.0.) The reason for this problem is due to the presence of 
small real eigenvalues in the fermion matrix. In an experiment, [9], we have 
chirally rotated the lattice quark mass, am q , away from the real axis; the 
same configuration then gave a well behaved pion propagator. (We might 
then expect some mixing in the correlation function of the particle with its 
parity partner, however for the pion in the quark model this partner does 
not exist.) So perhaps simply throwing away the configuration does not 
affect the spectrum (?). It is desirable to have a crude indicator of whether 
we have an exceptional configuration. In [10] the simple proposal was made 
to look at the pion norm, 

n({t/}) = E \j 5 G(x, t; 0, 0; {U}h 5 \ 2 . (3) 

x,t 

(In our application, the source (0, 0) was also Jacobi smeared.) In Fig. 4, we 
show a sequence of pion norms for j3 = 6.2. To decide on a criterion for an 
exceptional configuration (ie a spike in the pion norm) is not so easy. Some 
are obvious, for example from the pictures we have at k = 0.1354, n con f = 
217 a problem. (The corresponding pion propagator is shown in [9].) Closer 
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Figure 4- The pion norm, eg. (3), against configuration number for [3 — 6.2. 



to the critical point than here more spikes have been seen in Wilson data, 
[11]. To be safe, we have actually chosen a more conservative local criterion 
where if at any t value the (pion) correlation function fluctuates more than 5 
standard deviations from the local average we reject the configuration. This 
leads at = 6.2, 32 3 x64 for k = 0.1352, 0.1354, 0.13555 to rejection rates of 
about 2, 4 and 6% respectively. (For the latter two n values about 15% and 
33% of these rates were actually due to non-convergence of the inverter.) 
So in conclusion: while we feel that for any lighter quark mass than those 
considered here exceptional configurations become a real problem, at our 
masses while they are a nuisance, they do not distort the numerical result. 

We now turn to a consideration of the decay constant f ps . From eq. (1) 
we expect that af ps has no quenched chiral singularity in it, while a 2 gp 
diverges in the chiral limit. In Fig. 5 we plot the unrenormalised af ps . The 
results seem smooth over the whole quark mass range, with no singular 
behaviour. Looking at the ratio af ps /a 2 gp we see the same behaviour as in 
Fig. 1, ie for smaller quark masses than the strange quark mass a bend is 
seen in the data. As af ps /a 2 gp oc [(am ps ) 2 ] s then taking the logarithm gives 
a direct estimate of 5. In Fig. 6 we show this, with fit values 5 ~ 0.10(3) 
{(3 = 6.0) and 5 ~ 0.05(2) ((3 = 6.2) consistent with the previous results. 

Finally we consider chiral extrapolations of the nucleon and rho masses. 
In Figs. 7, 8 we show the results, together with a phenomenological fit 
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Figure 5. Unrenormalised (improved) af ps versus the pseudoscalar mass. 
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Figure 6. Unrenormalised (improved) ln[af ps /a 2 gp] against ln[((ro / a) 2 x (om ps ) 2 ]. The 
left line is a linear fit to the results for quarks with mass below the strange quark mass. 
The same notation as for Fig. 3. 



In distinction to eq. (2) this does not have a quenched linear chiral term. 
(As there is curvature in the results above the strange quark mass, we have 
considered (am p ^) 2 rather than am P: N and included a cubic term in eq. (4). 
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Figure 7. am p versus the pseudoscalar mass. The same notation as for Fig. 1. The fit 
function is given in eq. (4). 




Figure 8. amjv versus the pseudoscalar mass. The same notation as for Fig. 7. 



This gave a better fit function for the data.) While, for the nucleon this 
gives a good description of the data over the whole quark mass range and 
it is thus difficult to say in this case whether a linear term is necessary or 



9 



not, the p data might be showing some deviations for small quark masses. 
4. Conclusions 

Our main conclusion is that in quenched QCD there seems to be a danger- 
ous region for quark masses m q < m s . If we are interested in the strange 
quark mass or particles such as mjc, rnx*, ■ ■ -, or decay constants such as 
Ik, fx*, • • -j this does not represent a problem. For quantities involving only 
the u and d quarks, it is probably best to adopt the pragmatic approach of 
making fits for m q > m s and then to extrapolate this to m q ~ ?rj M /d, ie to 
chiral limit. (See, for example, the results for the quark mass, [7].) Evidence 
for chiral quenched logarithms is mixed - the best signal seems to be for 
the pion and its associated decay constant. Other channels seem to be less 
unambiguous. Indeed, as detecting and measuring quenching effects can be 
quite difficult this would indicate that the quenched approximation often 
seems to be working quite well. 

The above results should be regarded as preliminary. We hope to present 
full results shortly, [12], including continuum extrapolations (considered in 
the talk, but not described here). 
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